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ABSTRACT 

We have assembled a sample of halo stars in the solar neighborhood to look 
for halo substructure in velocity and angular momentum space. Our sample (231 
stars) includes red giants, RR Lyrae variable stars, and red horizontal branch 
stars within 2.5 kpc of the Sun with [Fe/H] less than -1.0. It was chosen to 
include stars with accurate distances, space velocities, and metallicities as well 
as well-quantified errors. With our data set, we confirm the existence of the 
streams found by Helmi and coworkers, which we refer to as the H99 streams. 
These streams have a double-peaked velocity distribution in the z direction (out 
of the Galactic plane). We use the results of modehng of the H99 streams by 
Helmi and collaborators to test how one might use velocity information and 
radial velocity information to detect kinematic substructure in the halo. We 
find that detecting the H99 streams with radial velocities alone would require a 
large sample (e.g., approximately 150 stars within 2 kpc of the Sun and within 
20° of the Galactic poles). In addition, we use the velocity distribution of the 
II99 streams to estimate their age. From our model of the progenitor of the 
H99 streams, we determine that it was accreted between 6 and 9 Gyr ago. The 
H99 streams have [a/Fe] abundances similar to other halo stars in the solar 
neighborhood, suggesting that the gas that formed these stars were enriched 
mostly by Type II supernovae. We have also discovered in angular momentum 
space two other possible substructures, which we refer to as the retrograde and 
prograde outliers. The retrograde outliers are likely to be halo substructure, but 
the prograde outliers are most likely part of the smooth halo. The retrograde 
outliers have significant structure in the direction and show a range of [a/Fe], 
with two having low [a/Fe] for their [Fe/H]. The fraction of substructure stars in 
our sample is between 5% and 7%. The methods presented in this paper can be 
used to exploit the kinematic information present in future large databases like 
RAVE, SDSSII/SEGUE, and Gaia. 
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Subject headings: methods: statistical, Galaxy: halo, Galaxy: kinematics and 
dynamics, solar neighborhood 



1. Introduction 

The kinematics of stars in the Galaxy provide information about its structure and for- 
mation history. Early Galaxy formation models, most notably those by Eggen et al. (1962), 
postulated a monoht hic collapse of gas an d dust, with a large cloud of gas collapsing and 



forming stars. Later, ISearle fc ZinnI (119781 ) suggested that the Galaxy was formed through 
a series of accretion events, with smaller proto-galaxies coming together to form a larger 
structure. Their scenario is more in keeping with our current u nderstanding of how struc 



tures formed in the ea rly universe in a "bottom-up" fashion (jSteinmetz fc Mulleii Il994 



White fc Springelll2nnnl ). 



Ibata et al. (1994) discovered that the Galaxy is currently accreting the Sagittarius 
dwarf spheroidal galaxy. Stars from this system are not yet well-mixed with the rest of 
the Galactic halo either spatially or in velocity space. Could other halo accretion events be 
detected if the stars from the accreted object have become wel l-mixed with the st ars in the 



halo spatially, but not in velocity or angular momentum space? iHelmi et al.l (1 19991 . hereafter 
H99), found evidence in angular momentum space for streams in the solar neighborhood 
using a sample of 97 metal- deficient red giants and RR Lyrae stars within 1 kpc of the Sun 
and, using the sample properties, generated a model of the streams. 

The purpose of this work is to search for further evidence for substructure in the local 
Galactic halo and to develop novel ways to detect halo substructure. In order to detect the 
subtle signs of kinematical substructure, accurate data are needed. We use two different halo 
star samples: the original H99 sample and a combined data set that includes local giants, 
red horizontal branch (RHB) stars, and RR Lyrae stars, all selected without kinematic bias. 
Both samples have full space velocities, distances, and metallicities. Further improvements 
to the solar neighborhood halo sample will be presented in the next paper in this series, 
H.L. Morrison et al. (2007, in preparation, hereafter M07), which will focus on the overall 
properties of the sample and their constraints on the formation history of the local halo, 
rather than specifically on substructure. 

Full space velocities are difficult to obtain for many stars, preventing a complete analysis 
of their motions like that of H99. However, the H99 models and original data show structure 
even in their radial velocity distributions, especially when stars in a particular direction, e.g. 
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near the Galactic poles, are considered. In addition, other possible substructures that we 
discuss in this paper occupy the tails of the angular momentum distribution, and so may 
be detectable in radial velocities dominated by rotational velocity. Standard statistical tests 
such as the Shapiro- Wilk test for normality, which is quite sensitive to the behaviour of the 
tails of the distribution, can be used to detect such deviations. We use the H99 model of the 
streams as well as a model of a smooth halo to show how streams can be detected in one 
component of velocity (f^), and we determine the detection limits for this method. We also 
extend this method to radial velocities for stars in the direction of the Galactic poles. 

In recent years, increasingly large and precise data sets of abundances for multiple el- 
ements have made it possible to study a new dimension of substructure in the halo: the 
chemical patterns. Abundances trace the star formation history of the objects that be- 
ca me today's halo. A good gen e ral di scussion of the possibilities of this technique is given 



m 



Freeman fc Bland-HawthornI (120021 ). In particular, the ratios of the a-elements such as 



Mg, Ca, and Ti to Fe allow the possibility of distinguishing between a chemical history 
where enrichment only comes from short-lived massive stars (Type II supernovae), and one 
where the contribution of Type la supernovae need to be considered as well. Type II su- 
pernovae preferentially enrich the interstellar medium (ISM) with a elements, while Type 
la supernovae produce more iron but need a longer period of chemical evolution. Thus the 
bulk of the local halo stars have a higher value of [a/Fe] than local disk stars, reflecting 
different formation conditions. However, there are some metal-poor s tars that have unusu- 



ally low values of [a/Fe], such as those in dwarf spheroidal galaxies (IShetrone et al.l 12001 



Kindll997 : 



20031 : iTolstov et al.ll20oi). a small nu mber of outer halo globular clus ters such as 



and Rup 106 (Brown et al. 1997 . 1999), and some local halo field stars JCarnev et al 



Nissen fc Schuste: 



1997 



Fulbrightll2002l : [Stephens fc Boesgaard 



2002 



'al 12 



1997 



Ivans et al. 



20031 : IVenn et al.ll2004l ). For stars from substructures identified via kinematic methods, we 
investigate whether there are any unusual patterns in [a/Fe]. 

Section [2] describes the compilation of the data sets used in this work in detail. Informa- 
tion about our model of the Galactic halo is given in Section [3] and our results are described 
and discussed in Section HI Section [5] summarizes our conclusions. 



2. Data Sets 

The search for substructure in H99 is particularly powerful because estimates of all three 
space velocities are available for their sample. This allows the use of angular momentum 
(approximately conserved in a roughly spherical potential) to isolate stars with similar ori- 
gins. In the pre- Gaza era, the need for three space velocities limits us to a relatively small 
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volume surrounding the Sun, where sufficiently accurate proper motions are available. It is 
important to note that methods of ide ntifying satellite debris using physical quantities such 



as an gular momentum and energy (e.g. iLynden-Bell fc Lynden-Belllll995l : iHelmi et al.lll999 



20061 ) rely strongly on good distance estimates. (The angular momentum estimate varies as 



distance squared, since two out of three velocity coordinates are obtained by multiplying the 
proper motion by the distance.) Thus the accuracy of distance estimates is very important. 



2.1. H99 Data 



The H 99 sample stars were se lected from the compilations of iBeers &: Sommer-Larsen 



( I1995I ) and IChiba fc Yoshiil (Il998l . hereafter CY98) to have distances less than 1 kpc and 
[Fe/H] less than -1.6. This sample contains nearby stars in three evolutionary states: ffist 
ascent red giants, RHB stars, and RR Lyrae variable stars. These stars present differ- 
ent challenges for stellar population work. Good velocity estimates are difficult to obtain 
for RR Lyrae stars, but are relatively straightforward for red giants. This is because RR 
Lyrae stars show significant velocity variation due to pulsation, while ffist ascent giants and 
RHB stars have spectra that lend themselves readily to accurate velocity estimation. Con- 
versely, distance estimates for RHB and RR Lyrae stars are good, but distance estimates 
for red giant stars less so. RR Lyrae and RHB stars show little variation in absolute mag- 
nitude with metallicity. For examp le. iMartin fc Morrison! (119981 ) estimated a typical error 
of 10% for RR Lyrae distances, and Vivas et al. ( 20051 ) estimate 6%. First ascent red giants 
generally have less precise distances, because the position of the giant branch in the color- 
magnitude diagram depends strongly on metallicity, and so metallicit y measurement errors 



propa gate to a larger distance error of order 20% (see discussions in [Morrison et al.lll990 



2OO3I ). The advantage of smaller distance errors for RHB stars is offset by the difficulty of 
identifying these stars; in order to separate the subtle differences in gravity between ffist 
ascent red giants below the horizontal branch and horizontal br anch stars we need pho- 



tometry from intermediate-band system s like Stromgren or DDO (lBondlll980l : iNorris et al. 



1985 : Anthonv-Twarog fc Twaroel 1994 ). Luckil y, many of the stars in the H99 sample were 
observed by lAnthony-Twarog fc TwarogI (11994 . hereafter ATT94) who obtained accurate 
Stromgren uvby colors, enabling them to classify stars as either RHB or ffist-ascent giants. 
Below we calculate the actual errors on distances, space velocities and angular momenta for 
the ATT stars, illustrating that these quantities are known remarkably precisely for a sample 
including red giants. 
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2.2. Combined Data Set 

Since H99 was publislied there have been several improvements on the infor mation avail- 



able fo r the solar neighborhood giant and horizontal branch stars. In particular, iBeers et al. 



( 120001 . BOO hereafter) presented a catalog of 2016 stars selected without kinematic bias, and 
we use preliminary results from a new sample of RR Lyraes with significantly improved 
radial velocity measurements (M07). 



2.2.1. RR Lyrae Variables 



While the iLaydenl (Il994l ) data set provided accurate [Fe/H] and distance measurements 
for the local RR Lyrae stars, the velocity measurements were less accurate (typical errors of 
30kms-^). MOT obtained more accurate velocities for the RR Lyrae stars accessible from the 
north, observing most stars more than once and using high-quality light curves obtained by 
T.D. Kinman (1991, private communication) to correct for the large radial velocity amplitude 
of eac h star. Velocities were corrected for pul sation using the synthetic v elocity curves of iLiu 
(1l99ll ). Typical velocity errors are 15 km s~^. iMartin fc Morrison! (119981 ) published improved 
values of proper motion for RR Lyrae stars by averaging Hipparcos proper motions with 



accurate ground-based determinations from the USNO ACT reference catalog (lUrban et al. 



19981 ). We have selected a sample of 96 stars with [F e /H] less than - 1.0 fr om these data sets. 
Metallicities for these stars are from lLavdenI (1994 ) : iLavden et al.l (119961 ). The distances of 
stars in the M07 sample (derived using the iLaydenl (Il994l ) period-luminosity relation) range 
from about 0.5 kpc (RR Lyrae itself) to 2.5 kpc. Our sample presently only includes RR 
Lyrae stars of type ab. 



2.2.2. Red giants and RHB Stars 

We have selected a subset of the red giant stars in the BOO catalog for our expanded 
sample. In order to preserve the distance accuracy so important to the calculation of angular 
momentum, we critically compared the distance estimates given in the BOO catalog with those 
given by earlier work such as ATT94. This uncovered some surprisingly large, systematic 
differences in distance estimates. 

Figure [T] shows that the BOO and ATT94 distances agree reasonably well for distances 
greater than 1 kpc. For distances less than 1 kpc, the ATT94 distances are on average a 
factor of 2 larger than the BOO distances, a remarkably large amount for local, well-studied 
stars. 
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Both papers u se B — V colors and the globular cluster giant branch parameterization of 



Norris et al.l (Il985l ) to estimate absolute magnitudes for the giant stars. BOO assumed that 



all giants are on the first-ascent giant branch, while ATT94 used the Stromgren cl index 
to separate horizontal-branch and asymptotic giant branch (AGB) stars from first-ascent 
giants. AGB star contamination is expected to be small because of the short AGB lifetimes, 
but horizontal-branch stars are more common. Blueward of {B — V)o — 0.9, we expect 
significant numbers of RHB stars to appear in the sample. Assuming that these stars were 
on the first ascent giant branch, as BOO did, will lead to significant underestimates of their 
distance. This error will become larger for the bluest stars. Figure [2] shows that this is 
indeed the case, with BOO distances being on average a factor of 2 smaller than the ATT94 
distances for {B — V)o = 0.6. 

Thus we have restricted our sample to giants from BOO with {B — V)q greater than 0.9, 
red giants and RHB stars from CY98 with distances determined by ATT94, and RR Lyrae 
stars from the M07 sample. We have also restricted the sample to stars with [Fe/H] less than 
-1.0 and distances less than 2.5 kpc, and eliminated stars with thick disk kinematics (see 
Section SSI for details). The {B — V)q criteria led us to accept 24 red giants and to reject 
85 red giants from BOO that had otherwise acceptable distances and metallicities. Fourteen 
of the accepted stars were later eliminated from the final sample because of their thick disk 
kinematics (see Section H73|) . For stars in our sample within 0.5 kpc of the Sun, our distances 
agree with Hipparcos distances (lESAl 119971 ) within the errors. The only exception is HD 
135449, whose distance in Tabled] is twice the distance given by Hipparcos. In future versions 
of this sample (M07), we plan on using the Hipparcos distances for nearby stars. Metallicities 
in BOO come from a large variety of sources, but the ATT94 data (which dominate our 
combined sample of red giants and RHB stars) are significantly more homogeneous, being 
based either on uvby photometry (89 stars, an estimated [Fe/H] error 0.16 dex) or accurate 
spectroscopic abundances from the literature (68 stars). 



2.2.3. Error Estimates 

Because the data in our combined sample come from a variety of sources and include 
stars in different evolutionary stages, it is preferable to calculate the errors on derived quan- 
tities individually, since they can vary significantly (see Table 1 and Fig. [TTl) . We prop- 
agated the known errors on distance, velocity, and proper motion through the calculation 
of VR,Vtf),Vz, Jz, and J± using a Monte Carlo calculation, drawing new values of the input 
quantities from Gaussian distributions, recalculating the derived values and measuring the 
standard deviation of their distributions. 
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For the RR Lyrae and RHB stars, distance errors are relatively small, and have been 
assumed to be 7%. For the first ascent red giants in the sample, it is more difficult to quantify 
errors because the metallicity is the main driver of distance errors and many authors do not 
quote metallicity error^. We have calculated distance errors individually for stars in the 
ATT94 sample (where accurate metallicity errors are given, and in fact a significant number 
of stars have multiple high-dispersion metallicity estimates) and for a small number of red 
giants, most of which are part of the outher groups discussed in Section 14.31 To calculate 



distance errors, we used B — V color and the globular cluster giant branch loci of iNorris et al. 



(119851 ) with a Monte Carlo calculation of the effect of metallicity and color errors on the 
estimate of absolute magnitude. 

For the five red giant stars not in the ATT94 sample, we made the following assump- 
tions. Four (HP 18 7 10, HP 174578, HP 214925 a r id CD -68 1881) have [Fe/H] values from 
Norris et al.l (119851 ). iTwarog fc Anthony- Twarod (119941 ) have shown that the PPO [Fe/H] 
calibration used for these stars has systematic errors in some metallicity ranges. We have 
assigned [Fe/H] errors of 0.25, 0.50, 0.20 and 0. 5 dex respectively to account for this. BP 
-|-30 2282 has metallicity from iHartkopf fc Yosd (1l982l ): we have assigned an [Fe/H] error of 
0.35 dex f or this star. There was also one giant in ATT94 (HP 128279) that was too blue 
to use the iNorris et al.l (119851 ) distance calibration for the Monte Carlo calculation. Luckily, 
it has a Hipparcos parallax of reasonable accuracy so we have simply used the parallax and 
its error rather than our other distance error estimation procedure. 



2.2.4- Final Sample 

We have error estimates for distances and derived quantities (velocities and angular 
momentum) for the great majority of the sample: 210 out of 231 stars. The combined 
sample has a median distance of 1.1 kpc, median distance error 7%, and median metallicity 
[Fe/H] = -1.7. Median errors on vr,v^,Vz velocities are 15, 20 and 11 km s'^^ respectively, 
making this sample well-suited for careful investigations of substructure. 

Table [1] lists values of distances, metallicities, radial velocities, and galactocentric ve- 
locities for the red giants and RHB stars in the combined sample as well as the associated 
errors for these quantities. The RR Lyrae data will be published in M07. In Table [H the 
radial velocities are heliocentric and the vr, v^, and velocities are relative to the cen- 



^Note that CY98 calculated distance errors simply by assuming that the error on each metallicity value 
in their sample was 0.16 dex, while ATT94's actual quoted errors vary from 0.01 to 0.31 dex. We have set 
errors quoted at less than 0.1 dex to be 0.1 dex in our calculations. 
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ter of the Galaxy on a left-handed coordinate system. This has R pointing away from the 
center, increasing in the direction of Galactic rotation, and z pointing towards the north 
Galactic pole (NGP). Note that our angular momenta are calculated on this left-handed 
system in order to directly compare our results with the results of H99. (See the Appendix 
for more details.) The velocities were corrected for the motion of the Sun and LSR using 
= 220.0 km s-\ Uq = -9.0 km s'K Vr^ = 12.0 km s-\ and Wc^ = 7.0 km s~^ on the 



Vl 



standard left-handed coordinate system (IBlaauw fc Schmidtlll965l : iMihalas fc Binneylll98ll ). 
The Sun was assumed to be 8.0 kpc from the Galactic center. See the Appendix for the 
transformation between the local. Sun-centered coordinate system and the galactocentric 
coordinate system. 

The sample was selected without kinematic bias, but is spatially incomplete: it has 
fewer stars with low \b\ because of the difficulty of identifying these rare halo stars in heavily 
reddened regions. We show the distribution of distance and \b\ in Figure [31 It can be 
seen that the RR Lyrae sample is more complete in |6| than the giant sample, although 
neither are fully complete in The giant sample has a lower mean distance than the RR 
Lyrae sample because of the intrinsic rarity of the RR Lyrae stars. Since all three velocity 
components can be measured for each star, the spatial incompleteness does not cause a bias 
in the velocity distribution directly. However, the bias toward high Galactic latitude does 
select agai nst stars with a fla t tened space distribut ion, such as the metal-poor stars in the 
thick disk (INorris et al.lll985l : Morrison et al.lll990l ). Therefore, we would expect any such 
stars to be under-represented in our sample. 

Since there will be a signiflcant num ber of binaries in an y sample of stars (for example, 
17% in the sample of metal-poor giants of lCarney et al.ll2003l ). we need to consider the effect 
of binarity on our estimates of distance and velocity. Fortunately, because our sample is 
composed of giants and horizont al branch stars , the e ffect of any undetected binaries on the 
distance estimates will be small. ICarney et al.l (120031 ) studied 91 metal-poor red giants and 
RHB stars, (68 in common with our sample) and identified 8 spectroscopic binaries. The 
median velocity amplitude of these binaries is 8 km s~^, which is small compared to our 
typical error on a U,V or W velocity (10-15 km s~^). While we expect a similar number 
of undetected spectroscopic binaries in our sample, the effect on our conclusions will be 
negligible. Similarly, the "velocity jitter" that is seen in some stars close to the giant branch 
tip is of an even smaller magnitude (of order 5 km s~^) and will not affect our conclusions 
either. 
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Model Information 



Helmi fc Whitd (Il999l ) described models used to study the disruption of satellite galaxies 
by the Milky Way. H99 used these to model the progenitor of the star streams they detected 
in the solar neighborhood. The remnants of a progenitor with an initial velocity dispersion 
of 18 km s~^ and core radius of 0.5-0.65 kpc fit the star streams discovered in H99. Its 
orbit has an apocenter at 16 kpc, pericenter at 7 kpc, a maximum height above the plane 
of 13 kpc, and a radial period of 0.4 Gyr. The present paper uses an improved version of 
the model, which better accounts for the self-gravity of the satellite and has positions and 
space velocities for 10^ particles in the disrupted progenitor from Gyr- 13 Gyr. After 3-5 
Gyr, these particles fill much of the space between 7 and 16 kpc from the Galactic center 
(the amount of debris from the progenitor peaks at the solar circle) and are very well-mixed 
spatially within ~5 kpc from the Sun. The disrupted progenitor still has some structure in 
velocity and angular momentum space, however, as can be seen in Figure |H 

In addition, we have created a model of a smooth halo. This model allows us to 
test the null hypothesis that the halo has no substructure. The density dist ribut i on of 
the smooth halo ( number of stars per unit volume) is proportional to r~^'^ ( IZinnI Il985 



Vivas fc Zinnll2006h. The ve locity distribution of the halo was modeled using a velocity el- 
lipsoid (jSchwarzschildl 119071 ). The models were created by randomly selecting points from 
these distributions. Figure shows one realization of this model in velocity and angular mo- 
mentum space. The parameters {vr), {v^), {vz), ctr, cr^, and were determined by finding 
the average and the standard deviation of the vr, v^, and Vz velocity distributions from the 
combined data set excluding known stream stars. 



The standard velocity ellipsoid numbers for the halo (e.g. iChiba fc Beerd |2000| ) were 
not used because the H99 star streams, in particular, bias the velocity distributions. For 
example. Figure [6] shows that the stream stars have large absolute z velocities. Including 
these stars in the velocity ellipsoid calculations changes o"^ from 84 km s~^ to 101 km s~^. 
Table [2] gives the derived velocity ellipsoid parameters for our combined sample (stars with 
distances less than 2.5 kpc and [Fe/H] < —1.0, excluding likely thick disk stars), both with 
and without stream stars. In the case of {vr) and (vz), the paramet ers used are consistent 
with zero within the errors. However, (v^) is not consistent with zero. IChiba &: Beerj (120001 ) 
find a similar (w^) for their sample within 1 kpc of the Sun and note that (v^) decreases 
when larger distance ranges are included. Note that excluding the star streams from the 
sample used to determine the velocity ellipsoid also increases the anisotropy of the halo. 
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4. Identifying Structure in the Tails of the Velocity Distribution 

4.1. Method 

Using distances and full space velocities of metal-poor red giants and RR Lyrae stars, 
H99 identified the remnant of a merger between a small satellite galaxy and the Milky 
Way. Since obtaining proper motions of stars is difficult, another way to identify subtle 
merger remnants is desirable. Here we explore signatures visible in the tails of the velocity 
distribution. Figure M plots the H99 data in velocity and angular momentum space. Its 
bottom panel plots Jz, the component of the angular momentum (per unit mass) out the 
plane of the Galaxy's disk, versus J_l = + Jy, the angular momentum (per unit mass) 
in the plane of the disk. The H99 stream stars {large circles) clearly occupy a different 
region of phase space than the rest of the halo stars, attesting to the debris' common origin 
as a Galactic satellite. Figure [6] {top panels) shows cylindrical velocity coordinates. The 
stream stars {large circles) are located in two clumps, one with Vz ~ —200 km s~^ and the 
other with Vz ~ 200 km s~^. The bimodality in Vz seen for the nearby stream stars comes 
about because their orbits reach to a large distance (13 kpc) above and below the Galactic 
plane. Therefore, any nearby substructure stars will be passing through with high positive 
or negative Vz values. We can anticipate that this bimodality will be less marked as more 
distant stars are included in the sample, because it will include stream stars with smaller 
values of \vz\. 

Note that the Vz components of the H99 streams extend well into the tails of the ve- 
locity distribution, as can be seen in Figure [71 A smooth halo velocity distribution is well- 



appro ximated as a multi-dimensional normal (Gaussian) distribution (e.g. lBinney &: Merrifield 



19981 ). so testing the Vz distribution for normality may reveal the presence of substructure 
like that in H99. In addition, the line of sight component of velocity for stars near the 
Galactic poles is dominated by the Vz component of the star's velocity. Therefore, if one 
uses a sample of stars near the Galactic poles, one should be able to test for structure in 
the Vz distribution using only radial velocities. While with this technique we are only able 
to identify the presence of substructure in the data, not necessarily the substructure itself, 
it will be useful in identifying data sets where other observers might want to aggressively 



pursue proper motions and distances^ 



^One might imagine employing a genetic- algorithm technique to find stream stars where random groups 
of stars are selected, the Shapiro- Wilk product calculated, and the groups that have the highest p- values are 
bred together to create the next generation of groups to test. This is, of course, a computationally intensive 
method. 
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There are ma ny tests for deviations from norraality ; one of the most powerful is the 
Shapiro- Wilk test Jshapiro fc Wilklll965l : [D^Agostinolll986h . This test is based on the concept 
of probabihty plots, which plot the cumulative distribution functiorj^l using a transformation 
of the vertical axis that makes normally distrib uted data fall along a straight line. Examples 
of probability plots are shown in Chapter 1 of iD'Agostino fc Stephens! ( 1l986l ). The slope of 
this line gives an estimate of a for a normal d istribution. The Shapiro- Wilk test compares this 
slope with the sample standard deviation (IStephensI Il986l ). This technique is particularly 
sensitive to de viations from normality in the tails of the distribution. Code for this test 
( Royston 1995 ) can be obtained at StatLib0 

To test this method, we applied the Shapiro- Wilk test to the H99 Vz data. Here the 
null hypothesis is that the data have a normal distribution, and the probability value 
(p-value for short) indicates how often data with a normal distribution would produce this 
data set. Thus the smaller the p-value, the less likely it is that the data come from a normal 
distribution. (Note, however, that a high p-value does not mean that the sample distribution 
is normal, just that it is consistent with a normal distribution.) The results of this test are 
given in Table [3l We list the tested data set and the associated p-value. The entire sample 
fails the Shapiro- Wilk test at the 1% level. When the H99 stream stars and one other star 
(HD 124358) at (J^ ~ —1500, J_l ~ 2300) are excluded, however, the sample passes the 
Shapiro- Wilk test. The additional excluded star, indicated by a large triangle, has a very 
retrograde orbit and is well away from the rest of the data in angular momentum space; it 
may belong to another stream (see Section These results demonstrate the effectiveness 
of this method for detecting streams like those in H99, where the stars occupy the tails of 
one component of the velocity distribution. 



4.2. Modeling Stream Detections with Distance 

Having established that we are able to detect substructure using the velocity dis- 
tribution, now we use the H99 model of the debris and our model of the smooth halo to 
see (1) how changing the sample distance limit affects our ability to detect streams and (2) 
how close stars need to be to the Galactic pole to detect substructure using radial velocity 
information only. 

To characterize stream detections, we produced random samples with sizes ranging from 



•^The cumulative distribution function gives the fraction of data points that are less than or equal to a 
value as a function of value: CDF{x) — P{X < x), where x is the value and X is the data. 

"'http: / /lib. Stat. cniu.edu/apstat/R94 
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20 to 4000 stars, and distance limits of 1, 2, 3, 4 and 5 kpc. In each case, we produced 10,000 
different realizations of the sample. While the stream stars are spatially well-mixed up to 5 
kpc from the Sun, the relative number of stream stars to smooth halo stars drops slowly as we 
move away from the Sun; see Figure [Hi The number of stream stars in each distance-limited 
sample was scaled using the models so that stream stars made up 10% of the sample for a 
distance limit of 1.0 kpc. The maximum distance limit in our simula tions was 5 kpc, since the 



veloc ity ellipsoid of the halo may change beyond this distance (e.g. IVedel fc Sommer-Larsen 



19901 ). For each sample, we tested the Vz distribution for normality using the Shapiro- Wilk 
test. We calculated the p-value (significance level) for each of the 10,000 realizations of the 
sample, and then calculated the average p-value for each sample. 

Figure M shows the average p-value as a function of sample size for various distance 
limits in our Vz samples. P- values below 5% fail the Shapiro- Wilk test for normality. The 
number of stars needed to detect the stream with a 1.0 kpc distance limit is less than the 
sample size of H99 sample with this distance limit (101 stars), so it is unsurprising that we 
were able to detect the H99 streams using this method (see Table |8]). As the distance limit 
increases, more stars are needed in the sample to detect the presence of the H99 streams. 
For example, with a distance limit of 5 kpc, approximately 600 stars are needed to detect the 
H99 streams, but only 75 stars are needed with a distance limit of 1.0 kpc. There are two 
causes here. First, the fraction of stream stars in a sample decreases slowly as the sample 
distance limit increases (see Fig. [H]). Second, the velocity distribution of the stream changes 
away from the solar neighborhood, and this results in fewer stream stars with extreme values 
of Vz- (Since the high-energy debris from the progenitor is bound to the Galaxy, stars nearest 
the plane will have the largest Vz velocities. As the distance from the plane increases, the 
Vz velocities of the stars will decrease.) The result of both effects means that increasing the 
distance limit on a sample does not necessarily increase the detection probability, if it does 
not have enough stars for the stream to be detected in the sample with the larger distance 
limits. 

To use radial velocities to identify kinematic structure in Vz, we need to identify how 
close the sample stars need to be to the Galactic poles to have enough Vz velocity information 
to detect structure. To examine this, we generated five different distance limited samples 
{D < 1.0, D < 2.0, D < 3.0, D < 4.0, and D < 5.0 kpc) and applied four different Galactic 
latitude limits to each sample (|6| > 40, \b\ > 50, \b\ > 60, \b\ > 70, and |6| > 80). The 
stream stars were assumed to be well mixed spatially so we used the same normalization 
for the stream stars as in the previous set of simulations. We calculated the radial velocity 
of each star in each sample and then used the same algorithm as above substituting the 
radial velocities for the Vz velocities to determine an average p-value for the radial velocity 
distribution of a sample. The results of these simulations are given in Figure [TOl 
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Comparing Fi gures [9] and [TOl one sees that for lower Galactic latitude limits, much larger 
sample sizes are needed to detect streams. As the Galactic latitude limits get further from the 
Galactic pole, the width of the smooth halo radial velocity distribution increases, drowning 
out the presence of stream stars in the wings of the distribution. As we sample closer to 
the Galactic poles, the sample size needed to detect the H99 streams using radial velocities 
approaches that needed to detect the streams in Vz velocity space. In short, using the radial 
velocity test to overcome a lack of proper motion data requires a significant increase in the 
sample size. 



4.3. Streams in the Combined Data Set 

We used the combined sample described in Section 12.21 to investigate further the appli- 
cability of the Shapiro- Wilk test to real data. We concentrate here on looking for outliers in 
the Jz versus Jj_ distribution. In future work (M07), we will investigate structure within the 
versus J\ distribution, i.e., w ithin the angular m omentum distribution of the disk and 



the halo themselves, such as in Navarro et al 



mm . 



While H99 used two different distance limits in their analysis (1 and 2.5 kpc) in order 
to limit the effect of larger errors in tangential velocity, our direct calculation of errors on 
the angular momenta make this unnecessary. II99 also limited their sample to stars with 
[Fe/H] < —1.6 to exclude thick disk stars. We chose simply to exclude such stars using their 
position on the angular momentum plot. To be considered a thick disk star, the star needed to 
have a Jz between 1500 and 2500 kpc km s~^ and a J± less than 600 kpc km s~^. These limits 



were chosen to match the angular momentum distribution of stars in the iNordstrom et al. 



(120041 ) sample with [Fe/H] greater than -1.0 and heliocentric distances less than 2.5 kpc. 
In addition, the star's needed to be near 220 km s~^ with all other velocity components 
near zero. Note that a star traveling at the speed (220 km s~^) and position (8.0 kpc) of 
the LSR has a Jz of 1760 kpc km s~^. The J± limit reflects the velocity dispersion of stars 
in the disk. These criteria excluded 24 stars, of which seven had [Fe/H] < —1.6. Note the 
symmetry about ~ for the resulting halo distribution (Fig. [TT]) . The dumpiness within 
the distribution of halo stars will be investigated in M07. 

We have identified the H99 stream stars and two other groups of outliers in the combined 
sample. The properties of the H99 stream stars and the outliers are given in Tables IH and 
[6l Figure [11] shows where these stars fall on the angular momentum and cylindrical velocity 
plot. The H99 streams are clumped in angular momentum and cylindrical velocity space, 
but not in Galactic coordinates. The first group of outliers, on the retrograde (left) side of 
the angular momentum plot, might also be tidal debris. It has extent and isolation similar 
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to the H99 streams. It also forms a kinematically distinct group in the cyhndrical velocity 
plots (at (vd, ) ^ —300 km s~^) as well as being completely made up of low metallicity stars 



(see Table IVenn et al.l (120041 ) have already noted the chemical homogeneity of the extreme 
retrograde stars in the solar neighborhood). It is not, however, as tightly clumped in velocity 
space as the H99 streams. The second group of outliers is located on the prograde (lower 
right) side of the angular momentum plot and has kinematics similar to disk stars, but with 
higher (("W^) ~ 300 km s~\ {vr) ~ km s~\ (f^) ~ km s"^). The prograde outlier 
group is also very metal poor (see Table ED . Because of the closeness to the disk and thick 
disk region of the angular momentum plot, it is possible th at this group is rela ted to the 



disk or was accreted into the disk by dynamical friction as in lAbadi et al.l (120031 ). 



We can estimate the probability of stars from a smooth halo populating the above 
regions in angular momentum space using the smooth halo model described in Section [31 
We randomly selected 231 stars within 2.5 kpc of the Sun from this model and counted 
the number of stars in each region of the angular momentum plot where we see outliers 
in the data. The size of the angular momentum region was chosen to be large enough to 
enclose all the outliers and their errors. We did not include AS Cnc in the H99 streams 
region because we are not sure if it is really part of this group. (See Section 14.3. 1[ ) We 
also randomly generated errors for the points from the model based on the distribution of 
errors in the combined data set. The errors were modeled as normal distributions with mean 
zero, a{Jz) = 186 kpc km s~\ and <j{J±) = 107 kpc km s~^. This process was repeated 
for 100,000 trials. Table [7] gives the box parameters, the probability that there would be 
as many or more stars as there are possible stream stars in the box if the halo was entirely 
smooth, and the average number of stars in the each region for a smooth halo. 

For the H99 stream, we found that there was a very low probability of having 11 smooth 
halo stars in that region of the angular momentum diagram. This result supports the con- 
clusions of H99. In the retrograde outlier region, there is a 7% chance of having six or 
more stars in that region of the angular momentum diagram. This set of outliers is possibly 
another stream. The simulations tell a different story for the prograde outliers. The prob- 
ability for having at least three stars in this region of the angular momentum diagram is 
70%. In addition, changing the criterion used to exclude thick disk stars could eliminate the 
prograde group or fill in the region between the halo distribution and the prograde group. 
The prograde outliers are likely smooth halo stars (perhaps with some contribution from the 
thick disk as well) rather than part of a stream. 
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4-3.1. H99 Stream Properties 



The H99 streams make up 5% of our total sample. For a subsample of the combined 
data set with the same parameters as the H99 sample ([Fe/H] < —1.0 and distances less than 
1.0 kpc), the fraction is 9%. We detect 11 of 12 stars that H99 detected and possibly add one 
more star to the streams (AS Cnc). The missing H99 star (HD 214925) we group with the 
retrograde outlier group rather than the H99 streams. In Figure [TTl AS Cnc is located very 
far from the rest of the H99 group in angular momentum space (at Jz ~ 4000 kpc km s~^ and 
J_L ~ 3600 kpc km s~^) , but agrees with the rest of the stars in velocity space. Therefore, 
we cannot say for certain whether AS Cnc is part of the H99 group without better data 
on its distance and space velocities. We include AS Cnc in our list of members of the H99 
streams for completeness, not because we are certain it is a member. 



Pilachowski et al.l (Il996l ) and iFulbrightl (120001 . |2002| ) gave high-dispersion abundance 



analyses for three of the H99 stream stars in Table 4: HD 128279, BD+30 2611 and HD 
175305. Their [Mg/Fe] and [a/Fe] abundances are similar to other halo stars in the solar 
neighborhood, suggesting that the gas that formed these stars was enriched mostly by Type 
II supernovae. 



Fiorentin et al.l (120051 ) has performed a similar analysis to ours using only the BOC) data 



and finds seven additional stars in the H99 streams. Two of the iFiorentin et al.l (120051 ) 
stars (HD 214161 and BPS CS 22189-0007) we excluded from our sample because they 
were classified by BOO as giants and had {B — V)n less than 0.9; t hey are likely to have 
underestimated distances (see Section 12.2.21) . IFiorentin et al.l (120051 ) also include RZ Cep 
in their list of members of the H99 streams. RZ Cep is a type c RR Lyrae and thus not 
included in the M07 sample of RR L yrae stars; it co u ld be an additional member of the 
stream. Finally, five of the stars that IPiorentin et aD J2005h detect (BPS CS 22948-0093, 
BPS CS 30339-0037, BPS CS 29513-0031, BPS CS 29504-0044, and BPS CS 22876-0040) are 
classified by BOO as turnoff stars. We did not include these stars in our sample because the 
distance estimates used in BOO for turnoff stars were based on UBV-photometry and thus not 
able to deal with the evolution of turnoff stars up to the subgiant branch. This ambiguity 
intro duces additional uncertainty into the distance estimates for these stars (ISchuster et al. 
20041 . especially their Fig. 8). 



We tested the Vz velocity data for normality using the Shapiro- Wilk test. The results of 
these tests are given in Table [Hi When the H99 stream stars are excluded from the sample, it 
tests positive for normality. Comparing our results to the simulations in Section 142] it is not 
surprising that we were able to detect the H99 streams in this sample using only Vz velocities. 
We did not test the radial velocity distribution of this sample for normality because there 
are not enough stars in the sample for this test to detect structure. 
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We can use the observed asymmetry in the Vz velocity distribution of the H99 streams 
along with our model of them to estimate the how long ago the progenitor was accreted. 
Initially, all stars are bound, and as time goes by, the stars are released leading to the 
formation of streams. These streams will phase-mix, and after many gigayears, this mixing 
will have progressed so that half the stars will cross the Galactic disk in each direction. In 
Figure [TTl we see that the observed fraction of stars in the stream with negative Vz is 0.72 
(8/11). To estimate the age of the H99 streams, we selected 11 model stars within 2.5 kpc 
of the Sun from our simulation and determined the fraction in the stream with negative 
Vz- We repeated this 1000 times to build up a probabihty distribution. We then compared 
the observed fraction of stars (8/11 = 0.72) to this probability distribution to determine 
an age (Fig. [T2|) . We find that the observed fraction is matched by the mean simulated 
fraction for a progenitor that was accreted between 6 and 9 Gyr ago. Although we cannot 
determine a definite age, we can rule out the possibility that the stream is either very young 
(3 Gyr) or very old (12 Gyr). Future observations can improve this situation. By increasing 
the observed total number of stars in the streams, the width of the peak of the probability 
distribution in Figure [12] decreases allowing a more accurate determination of the age of the 
accretion event. 



4-3.2. Retrograde Outlier Properties 

The retrograde outliers do not add significant structure to the Vz velocity distribution 
(the sample passes the Shapiro- Wilk test if the H99 streams, but not the retrograde outliers, 
are removed), but they do have significant structure in the direction, as can be seen in 
Figure [TTl To see if the Shapiro- Wilk test is able to pick out these deviations from normality, 
we ran various samples of data through the Shapiro- Wilk test. The results are shown in 
Table [9] Excluding the retrograde group of stars (Fig. [TT], inverted triangles) causes the 
sample to pass the Shapiro- Wilk test (test positive for normality). Indeed, all the samples 
that pass the Shapiro- Wilk test exclude this group of stars, which are very far from the 
rest of the velocities. Samples excluding only the prograde outliers do not test positive 
for normality. Evidently removing the prograde outliers from the sample just removes stars 
from the tail of the normal distribution of the halo and does not remove substructure. 

We note that this group of stars has an extreme ly low value o f J^,. The 38 globular 



clusters with proper motion estimates summarized by iDinescu et al.l (Il999l ) , most of which 
are within 10 kpc of the Galactic center, have Jz values ranging from -664 to 2307 kpc km s~^, 
while the mean value of Jz for our retrograde group is -2500 kpc km s~^. However, larger 
samples of nearby proper-motion selected stars have produced stars with even more extreme 
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negative vel ocities, for ex ample Kapteyn's star group, which has a mean Jz around -4000 
kpc km s~^ (lEggenlll996l ). 



Three of the stars from the retr o grade substructure also have hi gh-di spersion abundance 
information from iPilachowski et al.l (119961 ) , iMc William et al.l (119951 ) , and iGratton fc Sneden 
(1l988h : HD 6755, HD 200654 and CD-24 1782. Interestingly, these stars show a range 
of [a/Fe]: one has roughly normal [a/Fe] for the local halo while the other two have low 
a/Fe] for their [Fe/ H]. Note tha t our re trograde group corresponds roughly to the "extreme 



retrograde" class of I Venn et al. (l2004l ). and we find a similar behavior (somewhat lower 



[a/Fe] than normal) with our improved kinematical measures. 



4-3.3. Discussion 

Including both the H99 streams and the retrograde outliers gives a substructure fraction 
of 7% for the halo. The fraction of stream stars in the local halo is therefore at least 5% 



and possibly as high as 7% in this sample. iGouldl (120031 ) produced an estimate of the overall 



amount of substructure in the local halo (its "granularity" ) using an updated version on the 
NLTT proper motion sample. This method uses the fact that dominant streams in the solar 
neighborhood will not show spatial substructure over such a small volume, but will produce 
correlations between different components of velocity (non-zero cross terms in the velocity 
dispersion tensor). These will be reflected in the proper motion distrib ution. The fr action 



of stream stars in our sample is consistent with the upper limit set in iGouldl (120031 ). that 
one stream can comprise no more than 5% of the halo. 

Chiba and Beers (2000) noted that when they increase the sample size of halo stars to 
3 times that of II99 (728 stars), they only find nine stream stars total. There are several 
reasons why they may not have found more stream stars. First, the box that they use in their 
Figure 15 to identify the stream is too small. Our Figure H] shows the distribution of stars 
in the model within 2.5 kpc in angular momentum and velocity space. This distribution is 
much larger than the box that Chiba and Beers (2000) selects. They also do not remove stars 
with thick disk kinematics {Jz ~ 1750 kpc km~^ and ~ kpc km~^) from their sample, 
thus decreasing the percentage of II99 stream stars. Note, however, that the clumps in the 
V(j, — Vz velocity space are still apparent in their diagram, so it is still possible to pick out the 
H99 streams in their larger sample. This example illustrates the need to look not only at a 
angular momentum plot, but also a velocity plot to confirm membership in a moving group. 
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5. Conclusions 

We have assembled a sample of halo stars in the solar neighborhood to look for sub- 
structure in velocity and angular momentum space. Our sample of 231 stars includes red 
giants, RR Lyrae variable stars, and RHB stars within 2.5 kpc of the Sun with [Fe/H] less 
than -1.0. It was chosen to include stars with well-quantified errors and accurate distances, 
space velocities, and metallicities. Understanding the errors in the measured and derived 
quantities, especially distance, is crucial for this work since they may distort any underlying 
substructure. 

With our data set, we confirm the existence of the streams found by H99, which we refer 
to as the H99 streams. These streams have significant structure in their velocity distribution 
in the z direction (out of the Galactic plane). We use the results of H99 to test how one might 
use Vz velocity information and radial velocity information to detect kinematic substructure 
in the halo. We find that detecting the H99 streams with radial velocities alone would require 
a large sample (e.g., approximately 150 stars within 2 kpc of the Sun and within 20° of the 
Galactic poles). We also use the structure in the velocity distribution of the H99 streams to 
estimate the age of this group. From our model of the H99 progenitor, we determine that 
the H99 streams' progenitor was accreted between 6 and 9 Gyr ago. 

We have also discovered, in angular momentum space, two other possible substructures, 
which we refer to as the retrograde and prograde outliers. For the retrograde outliers, there 
is a low probability of that region of the angular momentum diagram being occupied by 
six or more smooth halo stars. Based on this evidence, the retrograde outliers are likely 
members of a stream. The prograde outliers, however, are most likely smooth halo stars 
(perhaps with some contribution from the thick disk as well) rather than part of a stream. 
The retrograde outhers display significant structure in the direction. Samples excluding 
the retrograde outliers pass our test for normality in the direction. The fraction of stars 
in our sample that are stream stars is between 5% and 7%. 

For H99 streams, the [Mg/Fe] and [a/Fe] abundances are similar to other halo stars 
in the solar neighborhood, suggesting that the gas that formed these stars was enriched 
mostly by Type II supernovae. The retrograde outliers show a range of [a/Fe]: one has 
roughly normal [a/Fe] for the local halo while the other two have low [a/Fe] for their [Fe/H]. 
Note that our ret rograde group corresponds roughly to the "extreme retrograde" class of 



Venn et al.l (120041 ). and we find a similar behavior (somewhat lower [a/Fe] than normal) 



with our improved kinematical measures. 

Although we are not the first to note that stellar debris from the disruption of a satellite 



would have a double-peaked distribution in galactocentric radial velocity (see also lMeza et al. 
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2005 



H99 



the methods de veloped in this paper add to the toolbox of kinematic methods (e.g. 



Helmi et al.| 120061) being developed to exploit future large databa ses such as RAVE 



fISteinmetz et al.l 120061 ). SDSSII/SEGUE flBeers et al.ll2004f ). and Gaia flTuron et al.ll2005l ) 
to detect kinematic substructure in our Galaxy's ha lo. These tools, in conjunction with 



studi es of spatial over-densities in the Milky Way (e.g. IWillman et al.ll2002l : iBelokurov et al. 



20061 ) . will provide crucial answers to the puzzle of how our Galaxy was formed. 
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A. Transforming Local Heliocentric Coordinates to Galactocentric Model 

Coordinates 

Given the wide variety of coordinate systems used in Galactic astronomy, we explicitly 
derive the transformation between the local. Sun-centered coordinate systems (x', y', z') used 
in this paper and the Galactocentric coordinate system used in the H99 models {x,y,z). 
Figure [13] illustrates the two coordinate systems. The local velocity coordinate system is a 
left-handed coordinate system with the x' axis pointing away from the Galactic center, the 
y' axis pointing in the direction of Galactic rotation, and the z' axis pointing in the direction 
of the NGP. Note that the L^/ angular momentum vector points in the direction opposite 
of the z' axis. In other words, the Galaxy rotates clockwise when you look down from the 
NGP. The model coordinate system used in H99 is also left-handed. 

The position of a star with respect to the Sun is given by 

x' = -rfcos(6) cos(/) (Al) 
y' = rfcos(6)sin(/) (A2) 
z' = dsm{b) (A3) 

where d is the distance to the star in kpc, b is the Galactic latitude, and / is the Galactic 
longitude. The Galactocentric model coordinates are then 

X = 8.0 kpc + x' (A4) 
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y 

z 



(A5) 
(A6) 



where 8.0 kpc is the distance between the Sun and the Galactic center. 



To transform the observed space velocities to the Local Standard of Rest (LSR) frame, 
we use the equations 



Vx' = U + Uq 

Vy' = V + Vq- 

Vz' = W + Wq 



(A7) 
(A8) 
(A9) 



where U, V, and W are the space velocities directed toward the Galactic anti-center, toward 
the direction of rotation, and toward the NGP. We used the va lues Vi^r = 220.0 km s" 



—9.0 km s 



12.0 km s ^, and Wq = 7.0 km s ^ (IBlaauw fc Schmidtl Il965 



Mihalas fc Binneylll98ll ) to correct for the motions of the Sun and the LSR. These velocities 
are the velocities in the model frame, i.e., the following 



(AlO) 
(All) 
(A12) 



The angular momentum components (per unit mass) are then the cross products 

Jx = yvz - VyZ (A13) 
Jy = zv,^ - v-^x (A14) 
Jz = xvy - v^y (A15) 

and J_L is i^J^ + Jy. Note that these angular momentum products are calculated on a left- 
handed system. While this does not make a difference in Jj_, a left-handed Jz points in the 
opposite direction as a right-handed Jz- In other words, the left-handed Jz of a star near 
the Sun is 220 km s~^ x 8 kpc = 1760 kpc km s~^, while the traditional right-handed Jz of 
a star near the Sun is —1760 kpc km s~^ 
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Table 1. Combined Sample Red Giants and RHB stars (see Section [2^2] for details and 

M07 for the RR Lyrae data). 
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CD 


-68 1881 


134.0 


1.19 


0.33 


-1.81 


-62.26 


28.0 


153.98 


19.6 


-173.00 


50.7 


BOO 


HD 


165195 


-0.2 


0.65 


0.03 


-2.42 


-151.91 


6.1 


14.94 


7.2 


-30.55 


2.3 


CY98 


HD 


166161 


68.4 


0.20 


0.15 


-1.27 


-146.25 


9.9 


72.32 


26.6 


9.01 


1.2 


CY98 


HD 


175305 


-181.0 


0.12 


0.19 


-1.54 


32.51 


17.0 


128.72 


10.7 


-222.72 


30.1 


CY98 


HD 


174578 


-4.0 


1.10 


0.36 


-1.69 


-1.42 


10.0 


350.17 


37.9 


94.00 


28.3 


BOO 


HD 


184266 


-348.4 


0.22 


0.07 


-1.87 


298.37 


1.0 


-50.17 


10.4 


-98.31 


14.5 


CY98 


HD 


184711 


101.6 


1.01 


0.04 


-2.30 


-74.43 


2.4 


-0.15 


9.3 


-104.23 


4.4 


CY98 


HD 


232078 


-391.0 


0.83 




-1.61 


137.49 




-136.65 




-75.00 




BOO 


HD 


186478 


30.9 


1.02 


0.06 


-2.45 


-192.92 


9.2 


-137.86 


25.3 


-67.16 


6.0 


CY98 


HD 


187111 


-181.0 


0.61 


0.18 


-1.95 


137.06 


2.9 


24.79 


21.5 


-99.15 


28.6 


CY98 


BD 


-18 5550 


-126.2 


0.74 


0.08 


-2.84 


36.96 


5.1 


-98.39 


22.6 


-96.32 


12.6 


CY98 


HD 


190287 


135.0 


0.14 


0.23 


-1.09 


-136.43 


3.5 


129.59 


26.3 


-63.15 


1.9 


CY98 


HD 


195636 


-257.8 


0.59 


0.07 


-2.82 


-36.60 


14.4 


-141.82 


17.0 


166.43 


4.5 


CY98 


BD 


-17 6036 


19.2 


1.27 


0.09 


-2.70 


-200.00 


17.0 


17.60 


22.3 


52.15 


11.6 


CY98 


BD 


-15 5781 


-76.4 


1.93 


0.08 


-2.47 


16.04 


7.5 


26.02 


18.7 


-61.81 


12.6 


CY98 


BD 


-14 5890 


117.5 


0.93 


0.19 


-2.01 


-267.26 


30.9 


-19.80 


60.3 


-130.64 


13.5 


CY98 


CD 


-37 14010 


-200.0 


1.68 




-2.55 


244.98 




-94.42 




-33.00 




BOO 


HD 


200654 


-48.0 


0.46 


0.06 


-2.79 


382.45 


21.2 


-344.74 


34.4 


-225.92 


15.9 


CY98 


BD 


-03 5215 


-293.7 


0.73 


0.07 


-1.72 


147.76 


4.3 


-30.82 


5.8 


96.55 


7.1 


CY98 


HD 


204543 


-98.2 


0.72 


0.11 


-1.69 


-32.91 


9.8 


47.15 


14.9 


-5.24 


7.5 


CY98 


HD 


205547 


47.0 


1.10 




-1.85 


-190.90 




46.41 




103.00 




BOO 


HD 


206739 


-57.8 


0.57 




-1.58 


78.66 




115.51 




-54.51 




CY98 


BD 


-09 5831 


14.5 


1.90 


0.13 


-1.87 


-68.55 


13.2 


126.24 


17.0 


-56.01 


9.2 


CY98 


HD 


235766 


-314.0 


0.92 




-2.35 


-315.17 




6.82 




3.00 




BOO 


HD 


214362 


-92.3 


0.49 


0.07 


-2.20 


321.17 


21.0 


-20.30 


15.4 


-133.39 


15.4 


CY98 


HD 


214925 


-328.0 


2.15 


0.08 


-2.14 


93.56 


9.5 


-289.91 


36.1 


143.00 


14.0 


BOO 


HD 


216143 


-115.9 


0.69 


0.07 


-2.20 


-327.79 


25.1 


-29.90 


15.9 


76.35 


2.4 


CY98 


HD 


218857 


-169.6 


0.41 


0.19 


-2.15 


-129.95 


32.0 


44.21 


23.7 


153.77 


1.2 


CY98 


HD 


220662 


-78.1 


1.87 


0.08 


-1.59 


0.33 


11.2 


10.05 


17.2 


37.78 


5.0 


CY98 


HD 


220838 


-22.8 


1.43 


0.07 


-1.72 


-24.64 


8.4 


149.05 


8.0 


22.65 


3.1 


CY98 


HD 


221170 


-121.9 


0.69 


0.07 


-2.01 


-151.22 


8.6 


93.53 


3.3 


-62.75 


8.8 


CY98 
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Table 1 — Continued 



Name 


RV 


D 


S D 


[Fe/H] 


Vr 


5 Vr 


V4, 


5 v^ 




<5 Vz 


Source 




km s^^ 


kpc 


kpc 




km s^-'^ 


km s^^ 


km s^^ 


km s^-"^ 


km s^-"^ 


km s"'^ 




HD 222434 


13.4 


0.97 


0.14 


-1.56 


-98.84 


12.2 


43.08 


26.1 


10.87 


3.6 


CY98 
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Table 2. Smooth Halo Model Parameters 





With H99 Streams 


Without H99 Streams 


{vr) 


-0.93 ± 10.20 


-3.57 ± 10.63 


(JR 


154.96 ±7.21 


157.30 ±7.52 


{V4,) 


28.62 ±7.19 


23.42 ± 7.41 




109.33 ±5.09 


109.64 ±5.24 




-6.43 ±6.66 


-1.32 ±5.66 




101.23 ±4.71 


83.75 ± 4.00 



Table 3. Shapiro- Wilk Results for H99 data 



Data Set Ngtars p-value 



All Stars 

Excluding H99 Streams 

Excluding H99 Streams & HD 124358 



101 

94 
93 



0.009 

0.05 

0.19 
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Table 4. H99 Stream Stars in Combined Data Set 



Name 


RV 


D 5 D [Fe/H] 


Vr 




SVr 








5v^ 








SVz 




Source 




km 


kpc kpc 




km s 


-1 


km s~ 


■ 1 


km s 


-1 


km s~ 


■ 1 


km s 


-1 


km s~ 


■ 1 




CD -36 1052 


306.0 


0.72 0.07 


-2.19 


87. 


,54 


3, 


,9 


137, 


,36 


3, 


,7 


-267, 


,16 


4, 


,8 


CY98 


AS CNC 


258.0 


2.31 0.16 


-1.89 


16, 


,09 


39, 


,7 


38, 


,16 


54, 


,2 


363, 


,00 


49, 


,4 


M06 


TT CNC 


47.0 


1.22 0.09 


-1.58 


106, 


,04 


12, 


,0 


87, 


,25 


13, 


,8 


-231, 


,10 


22, 


,0 


M06 


TT LYN 


-67.0 


0.65 0.05 


-1.76 


124, 


,42 


13, 


,0 


106, 


,61 


8, 


,8 


-233, 


,30 


13, 


,9 


M06 


HD 237846 


-302.8 


0.82 0.10 


-2.67 


-199, 


,16 


3, 


.1 


116, 


,48 


4, 


,7 


-195, 


,83 


3, 


,1 


CY98 


BD +10 2495 


252.9 


0.78 0.17 


-2.14 


136. 


,83 


35, 


,0 


175, 


.60 


4, 


,0 


288, 


.68 


7, 


,3 


CY98 


HD 119516 


-287.0 


0.51 0.07 


-2.49 


145. 


,36 


5, 


.4 


143, 


.15 


7, 


,5 


-250, 


.47 


1, 


,6 


CY98 


HD 128279 


-81.7 


0.16 0.22 


-2.22 


-21, 


,13 


15, 


,4 


140, 


,32 


27, 


,3 


-252, 


,33 


49, 


,9 


CY98 


BD +30 2611 


-278.2 


1.11 0.10 


-1.26 


13, 


,91 


9, 


,0 


155, 


,60 


4, 


,6 


-276, 


,18 


4, 


,7 


CY98 


AR SER 


127.0 


1.71 0.12 


-1.78 


98, 


,02 


21, 


,3 


100 


,60 


24, 


,8 


313, 


,40 


23, 


,1 


M06 


HD 175305 


-181.0 


0.12 0.19 


-1.54 


32, 


,51 


17, 


,0 


128, 


,72 


10, 


,7 


-222, 


,72 


30, 


,1 


CY98 


XZ CYG 


-146.0 


0.50 0.04 


-1.52 


26, 


,64 


4, 


.5 


152, 


.18 


11, 


,9 


-232, 


,50 


15, 


,2 


M06 



Table 5. Retrograde Outlier Stars in Combined Data Set 



Name 


RV D (5 D [Fe/H] 


Vr 


5vr 


V4, 








Source 




km s ^ kpc kpc 




km 


km s^^ 


km s^^ 


km s^^ 


km s^^ 


km s^^ 




HD 6755 


-328.7 0.17 0.18 


-1.62 


207.21 


72.1 


-326.56 


50.3 


110.96 


16.7 


CY98 


CD -24 1782 


118.5 0.73 0.08 


-2.37 


-137.14 


16.8 


-276.76 


38.1 


8.50 


9.0 


CY98 


HD 124358 


325.0 1.13 0.12 


-1.98 


104.03 


38.7 


-290.72 


55.5 


313.68 


8.8 


CY98 


RV CAP 


-84.0 1.06 0.07 


-1.72 


-69.03 


9.5 


-279.92 


34.3 


-193.70 


18.3 


M06 


HD 200654 


-48.0 0.46 0.06 


-2.79 


382.45 


21.2 


-344.74 


34.4 


-225.92 


15.9 


CY98 


HD 214925 


-328.0 2.15 0.08 


-2.14 


93.56 


9.5 


-289.91 


36.1 


143.00 


14.0 


BOO 



Table 6. Prograde Outlier Stars in Combined Data Set 



Name 


RV D <5 D [Fe/H] 


Vr 


SVr 




5v^ 


Vz 


SVz 


Source 




km kpc kpc 


km s^-"^ 


km s^^ 


km s^-"^ 


km s^-"- 


km s^-"- 


km s^-"^ 




TV LEO 


-97.0 1.84 0.13 -1.97 


-113.72 


30.4 


326.94 


24.2 


-22.90 


20.9 


M06 


BD +01 3070 


-329.4 0.50 0.23 -1.85 


343.77 


27.3 


273.99 


12.6 


-112.58 


26.0 


CY98 


HD 174578 


-4.0 1.10 0.36 -1.69 


-1.42 


10.0 


350.17 


37.9 


94.00 


28.3 


BOO 
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Table 7. Probability of Finding Halo Stars in Various Regions of the Angular Momentum 

Diagram 







Lower Left 


Upper Right 






Group 


Number of Stars 


Corner 


Corner 


P{n > n*) 


<n> 


H99 Streams 


11 


[500, 1400] 


[1500,2500] 


0.00024 


2.81 


Retrograde Outliers 


6 


[-3000, 0] 


[-1750, 2500] 


0.069 


2.86 


Prograde Outliers 


3 


[2000, 30] 


[3000, 1200] 


0.69 


3.55 



Table 8. Shapiro- Wilk test results for combined sample Vz data {D <2.5 kpc and 

[Fe/H] < -1.0) 



N 


Prograde 
▲ 


Excluded Groups 
Retrograde H99 streams 
▼ • 


p- value 


231 








0.002 


219 






X 


0.584 


216 


X 




X 


0.584 


210 


X 


X 


X 


0.877 


213 




X 


X 


0.850 
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Table 9. Shapiro- Wilk test results for combined sample data {D < 2.5 kpc and 

[Fe/H] < -1.0) 



N 


Prograde 
▲ 


Excluded Groups 
Retrograde H99 streams 
▼ • 


p- value 


231 






0.00014 


228 


X 




1.9e-06 


222 


X 


X 


0.01 


225 




X 


0.21 


219 




X 


0.00061 


216 


X 


X 


l.le-05 


210 


X 


X X 


0.068 


213 




X X 


0.41 
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Fig. 1. — Difference between distance estimates for red giants from ATT94 and BOO as a 
function of distance. The differences are expressed as a fraction of the ATT94 distance. It 
can be seen that for stars closer than 1 kpc, the BOO distances are on average 40% smaller 
than the ATT94 distances. 
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Fig. 2. — Difference between red giant distance estimates from ATT94 and BOO, shown as a 
function of (B-V)o color. It can now be seen that the problematic stars are the bluest ones, 
with (B-V)o < 0.9. 
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Fig. 3. — Distribution of distance and absolute value of Galactic latitude (6) for (a) the RR 
Lyrae variables in our sample and (b) the red giant and RHB stars in our sample. 
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Fig. 4. — Structure in velocity and angular momentum space for the H99 star streams model 
for stars within 2.5 kpc of the Sun. The bottom panel is a plot of the angular momentum 
in the plane of the Galaxy's disk as a function of angular momentum out of the disk. The 
top panels show the cylindrical velocity coordinates of the data plotted against each other. 
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Fig. 5. — Structure in velocity and angular momentum space for the smooth model for stars 
within 2.5 kpc of the Sun. 
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Fig. 6. — Plots of the distribution of the H99 data in velocity and angular momentum space 
for stars with [Fe/H] < —1.6 and distances less than 1.0 kpc from the Sun. The stream stars 
are represented by the large circles. The additional Vz outher (HD 124358) is shown by the 
upward pointing triangle. 
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Fig. 7. — Histogram of the H99 data velocity distribution with the stream stars shaded. 
A normal distribution is plotted as a dashed hne. Note how the H99 stream stars widen the 
wings of the velocity distribution making it deviate from normality. 
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Fig. 8. — Fraction of stream stars in halo as a function of distance limit. Note that the 
fraction of stream stars decreases as the distance limit increases when the fraction of stream 
stars is normalized to 10% of the halo stars within 1 kpc of the Sun. 
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Fig. 9. — Average p-value for 10,000 random samples of Vz velocities with the relative number 
of stream to smooth halo stars fixed at 10% for a distance limit of 1.0 kpc. The standard 
deviation of the trial p- values is on the order of the width of the lines. The distance limits 
used for the samples are 1.0, 2.0, 3.0, 4.0, and 5.0 kpc. As the distance limit increases, the 
size of the sample for which detections are possible also increases. 
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Fig. 10. — Average p-value for 10,000 random samples of radial velocities (corrected to the 
Galactocentric frame of reference) with 10% stream stars for a distance limit of 1.0 kpc. The 
standard deviation of the trial p-values is on the order of the width of the lines. Galactic 
latitude limits > 40, |6| > 50, |6| > 60, |6| > 70, and |6| > 80) were imposed on samples 
with five different distance hmits (1.0, 2.0, 3.0, 4.0, and 5.0 kpc). N.B. The x-axes on these 
graphs all have different scales, since the range of samples tested is very different in each 
case. 
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Fig. 11. — Angular momentum and cylindrical velocity plots for the combined data sample 
with the selection criteria: D < 2.5 kpc and [Fe/H] < —1.0. The retrograde outliers are 
indicated by downwardOfacing triangles, the prograde outliers by upward-facing triangles, 
and the H99 streams by large circles. Note that the empty region of the angular momentum 
diagram at Jz = 1500 — 2500 kpc km s~^, J_l < 600 kpc km s~^ was occupied by likely 
thick-disk stars. These stars were excluded from our final sample based on their position in 
the angular momentum diagram. 
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Fig. 12. — Probability that the stream with negative Vz will make up a given fraction of 
the total number of stars from the accreted satellite near the Sun calculated using the H99 
model at four different times: 3, 6, 9, and 12 Gyr. The dashed line indicates the observed 
fraction of stars in this stream (8/11 = 0.72). 
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Fig. 13. — H99 model {x,y,z) and the local, Sun-centered {x\y',z') coordinate systems. 
Note that both the H99 model coordinate system and the local, Sun-centered coordinate 
system are left-handed. The direction of Galactic rotation is in the positive y'-direction. 



